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Abstract 

The two-fluid model of liquid helium is generalized to the case that the super- 
fluid fraction has a small entropy content. We present theoretical arguments 
in favour of such a small superfluid entropy. In the generalized two-fluid 
model various sound modes of He II are investigated. In a superleak carrying 
a persistent current the superfluid entropy leads to a new sound mode which 
we call sixth sound. The relation between the sixth sound and the superfluid 
entropy is discussed in detail. 
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1 Introduction 



The two-fluid model can be considered as the fundamental theory for the 
hydrodynamics of liquid helium below the A-point. One of the model as- 
sumptions is that the entropy of the superfluid fraction vanishes. A small 
superfluid entropy S s (of the order of one percent of the total entropy S) 
is, however, not excluded by the experiment. In this paper we generalize the 
two-fluid model to the case of a small superfluid entropy and investigate the 
consequences. 

The paper Entropy of the Superfluid Component of Helium by Glick and 
Werntz [|IJ finds that the ratio S s /S is less than 3%; earlier experiments 
cited in Ref. [l| are less accurate. A newer experiment which is sensitive 
to this point has been performed by Singsaas and Ahlers @. In this work, 
second sound measurements are interpreted as that of entropy. As a result, no 
difference to the true (caloric) entropy has been found. Due to the uncertainty 
in the absolute values of the caloric entropy (about 1 to 2%) the limit for 
S s /S is not much lower than that of Ref. 0. 

In nearly all theoretical approaches S s is taken to be zero. As discussed by 
Putterman H this is an assumption. Microscopically, S s = follows from the 
identification of the superfluid density p s with the square of a macroscopic 
wave function [[|. This identification is plausible but unproven. In section 1.2 
we present theoretical arguments in favour of a small superfluid entropy. 

The remainder of this paper is based on the hypothesis of a non-vanishing 
superfluid entropy of a size which is not excluded by the experiment. We 
generalize the two-fluid model for this case (section 2). In the generalized 
two- fluid model we determine various sound modes (section 3). We find a 
new sound mode (sixth sound) which exists only if the superfluid entropy 
does not vanish. Section 4 discusses in detail an experiment by which the 
sixth sound could be detected. If the 6th sound exists this experiment would 
determine the superfluid entropy S s . If the 6th sound does not exist this 
experiment would yield a considerably lower upper limit for S s . 

1.1 Two-fluid model 

We start with a short review of the two-fluid model. Based on ideas by 
Tisza the two-fluid model was developed by Landau || . We use Putterman's 
monograph || as a standard reference for this model. Without dissipative 
terms the two-fluid equations read 



d t p + V(p n u n + p s u s ) = , 



(1) 



dt(ps) + V(psu n ) = 0, 
d t (p n u n + p s U s )i + dj (PSij + p n u ni u nj + p s u si u sj ) = , 



(2) 
(3) 
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md t u s + m(u s V)u s = — V/i. (4) 

We use the abbreviations d t = d/dt and di = d/dxi. In eq. (3) we sum over 
the index j. The entropy per particle is denoted by s = S/N, the mass of a 
helium atom by m, the mass density by p = mN/V, the velocity field by u, 
the pressure by P and the chemical potential by fi. The indices n and s refer 
to the normal and superfluid phase, respectively. 

Eq. (1) represents the mass conservation, eq. (2) the entropy conserva- 
tion (for reversible processes), eq. (3) the momentum conservation or Euler 
equation, and eq. (4) the motion of the superfluid. The eight equations (1) 
- (4) describe the dynamics of eight independent macroscopic variables. As 
independent variables one might chose the temperature T, the pressure P 
and the velocities u n and u s . These variables are fields depending on the 
coordinate r and on the time t. 

All macroscopic quantities X are functions of the eight variables. Galilean 
invariant macroscopic quantities can be written as a function of three vari- 
ables only, 

X = X(T, P, w 2 ) where w = u n — u s . (5) 

Most results (like the sound velocities) are eventually expressed by equilib- 
rium quantities X(T, P, 0) depending on T and P only. 

The eqs. (1) - (4) are supplemented by dissipative terms and by Onsager's 
quantization rule. The resulting two-fluid model is the fundamental theory 
for the hydrodynamics of He II. We restrict ourselves mainly to eqs. (1) - (4) 
for which we discuss the modification due to S s ^ 0. 

The two-fluid model is a phenomenological macroscopic theory @. It 
is based on macroscopic conservation laws (like mass conservation and the 
first and second law of thermodynamics) and on experimental evidence (like 
S s ~ 0). There are, however, also theoretical ideas on a microscopic level 
which support these equations. We refer in particular to London's postulate 
0. U that the superfluid consists of a macroscopic number of particles mov- 
ing coherently in a single quantum state. This conception explains that the 
superfluid entropy vanishes which is implicitly assumed in (2). Furthermore, 
it implies that the superfluid velocity is a gradient field, or that 

curl u s = . (6) 

This statement is contained in (4). The condition (6) may be violated in 
a vortex. For a vortex line London's macroscopic wave function leads to 
Onsager's quantization rule. 

1.2 Superfluid entropy 

In this subsection we present arguments in favour of a small superfluid en- 
tropy S$. These arguments lead to a theoretical estimate of S s . 
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The non-existence of the A-transition in 3 He proves that this transition 
is due to the exchange symmetry of the 4 He-atoms. The interatomic forces 
are about the same in liquid 3 He and 4 He; they are not the cause of the A- 
transition. These facts suggest a close connection between the Bose-Einstein- 
condensation of the ideal Bose gas (IBG) and the A-transition of liquid 4 He. 
This point of view has been put forward by London and has subsequently 
been supported by several authors || |3|. 

The condensate of the IBG forms a macroscopic wave function and pro- 
vides thus a model for the superfluid motion. At this point there is, however, 
a serious discrepancy between liquid helium and the IBG. The critical be- 
haviour of the condensate fraction of the IBG is 

^oc|t^, = 1, (7) 
P 2 

where t = (T — T\)/T\ is the relative temperature. In contrast to (7) the 
superfluid fraction of He II behaves roughly like 

^oc|t|^, vh\. (8) 
P 3 

Just below the transition this implies po p s - The critical exponent (3 = 1/2 
is an essential feature of the IBG. It does not appear possible to change this 
value by some modification of the IBG. (Note that the IBG free energy is 
not a logical starting point for a renormalization procedure because it already 
applies to an infinite system.) 

We present now a possible scheme which reconciles (7) with (8) preserving 
at the same time the essential features of the IBG (like (3 = 1/2). 

Following Chester || we use the IBG wave function together with a Jas- 
trow factor F = Hfij] such an ansatz is based on Feynman's discussion ||. 
Allowing for a condensate motion the many-body wave function reads 

* = 5F[exp(i$)P II t^P- ( 9 ) 

Here S denotes the symmetrization operator. The (p^ are the real single 
particle functions of the non-condensed particles and the i%k are the oc- 
cupation numbers. The schematic notation [y?k] nfe stands for the product 
V 9 k( r i) • V 9 k( r 2) ■ • • • ■ V ? k( r n Al ); this notation applies also to [exp(i$)] n °. All n 
particles adopt the same phase factor exp(i$(r)) forming the macroscopic 
wave funtion 

^(r) = ^exp(i$(r)). (10) 

This implies that the condensate particles move coherently with the (small) 
velocity u s = V$/m. 
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Eqs. (9) and (10) are a standard description |3j for a superfluid motion in 
an IBG-like model. (Actually, one has to construct a suitable coherent state 
lOfl instead of (9). This point is, however, not essential for the following 



discussion.) In this description the superfluid fraction p s / p equals the con- 
densate fraction n /N = po/p. We note that the current po u s is not depleted 
by the real Jastrow factors (in contrast to the condensate density po itself). 

In order to dissolve the discrepancy between (7) and (8) we assume that 
non- condensed particles move coherently with the condensate. This is possible 
if non-condensed particles adopt the macroscopic phase of the condensate: 

M> = SF[exp(i$)p J] [^ k exp(i$)P [] [ ¥*]"*• ( n ) 

For the single particle states with nk 3> 1 this phase ordering requires only 
a very small entropy change. The ansatz (11) assumes therefore this phase 
ordering for low momentum states. The limit k c up to which the particles 
move coherently may be considered as a model parameter. 

We evaluate the quantum mechanical expectation value (\I/|j|\E') of the 
current operator j for (11). The result is proportional to u s = l«V$/m; nei- 
ther the real Jastrow factor nor the real single particle functions contribute. 
Equating the statistical expectation value of this current with p s u s yields 

p p N *tt 

For (nk) we use the occupation numbers of the IBG-form. Fitting (12) to the 
experimental superfluid density determines k c (t). In this way the discrepancy 
between (7) and (8) is removed. The asymptotic behaviour p s oc |t| 2//3 implies 
k c oc |t| 2 / 3 . 

With (12) fitted to the experimental superfluid density, the entropy of the 
contibuting non-condensed particles can be calculated. The resulting predic- 
tion Jll] of this superfluid entropy S s is shown in Fig. 1. Due to k c oc \t\ 2 / 3 
the superfluid entropy per particle, s s = S s /N s oc k 2 oc \t\^ 3 , vanishes for 
T — > T\, 

s s (T x ,P) = 0. (13) 

Because of po/p — > 1 it vanishes also for T — ► 0. 

The experimental superfluid fraction p s /p obeys rather well the law of 
corresponding states ||12|| ; this means that it can be written as a function of 
t = t(T, P) = T/T\(P) — 1 alone. Since s s is determined by a fit to p s /p this 
should also hold for the superfluid entropy, 

s s (T,P)^g(t)=g(T/T x (P)-l). (14) 
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Figure 1: Prediction | l]] of the superfluid density S S (T, -Psvp) as a function 
of the temperature and at saturated vapour pressure. 



In this way the full T- and P-dependence follows from the prediction shown 
in Fig. 1 (using the lower scale only). The superfluid entropy (14) is of the 
form (5) with w — 0. 

The presented argument for a small superfluid entropy is based on the ob- 
vious relation between the Bose-Einstein-condensation and the A-transition, 
together with the discrepancy between (7) and (8). This argument serves as 
a motivation of our investigation. The remainder of this paper could just as 
well be based on the mere hypothesis of a non-vanishing superfluid entropy. 

2 Modified two-fluid model 

In this section we generalize the two-fluid model to the case of a small su- 
perfluid entropy s s = s s (T, P,w 2 ). 

2.1 Mass, momentum and entropy conservation 

In eqs. (1) and (3) the normal and the superfluid part are treated in a sym- 
metric way. Therefore these equations are unchanged by a non-vanishing 
superfluid entropy: 

d t p + V(p n u n + p s u s ) = for s s t^O, (15) 

d t (p a u n + p s u s )i + dj(PSij + p a u ni u nj + p s u si u sj ) = for s s t^O. (16) 

In eq. (2) the entropy density ps is carried by the normal fraction alone; 
consequently, the entropy current density is psu Q . For s s ^ the entropy 
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density is partly carried by the superfluid fraction leading to a different en- 
tropy current density: 

psu Q — > p n s n u n + p s s s u s . (17) 

Here s s = S s /N s and s n = (S — S s ) / (N — N s ) . The entropy continuity equation 
reads now 

d t (ps) + V(p n s n u n + p s s s u s ) = for s s 7^0. (18) 



2.2 Equation of superfluid motion 

In the static limit eq. (4) becomes V/i = — s VT + v VP = where v 
V/N = m/p. This yields the well-known fountain pressure (FP) 



dT J 



- for s s = . (19) 



FP 



V 



The origin of the FP may be explained as follows: Consider two containers 
with liquid He II connected by a superleak; initially both containers have 
the same temperature T and pressure P. By increasing the pressure in one 
container some superfluid liquid is pushed through the superleak to the other 
container. Since the superfluid fraction carries no entropy the other container 
becomes colder; effectively dP > in the first container is accompanied by 
dT > 0. According to (19) the ratio dP/dT is proportional to the missing 
entropy s per transferred particle. A possible superfluid entropy s s = S s /N s 
diminishes the missing entropy per transferred particle, which means 

(entropy deficit of a superfluid , . 

particle relative to average). 

As in (17), this replacement reflects the change in the entropy transport due 
to s s ^ 0. The replacement (20) applies, however, to the static limit and does 
not immediately yield the wanted generalization of (4). 

A basic property of the superfluid motion is curlu s = 0, eq. (6). This 
property is well-established experimentally. Theoretically, it follows from the 
conception that the supervelocity is the gradient of a macroscopic phase; 
this conception is not altered by the modified picture presented in section 
1.2. Eq. (6) implies that the l.h.s. of (4) is a gradient field. Therefore, the 
generalization of (4) must be of the form 

md t u s + m(u s V)u s = -V(/x - p s ) for s s ^ . (21) 

The chemical potential p yields the contribution —dp/dT = s in (19). Con- 
sequently, the replacement (20) implies 

dpJT,P,w 2 ) . 
- V dT ' = Ss(T, P, w 2 ) . (22) 
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Together with s s , eq. (13), the discussed modifications should vanish at T A . 
Therefore 

Ms (T, P, w 2 ) = - [ T dT s s (T', P, w 2 ) . (23) 

Eq. (21) with (23) defines the generalization of (4). 

We determine the FP from the generalized equation of motion (21). The 
FP experiment is done for u n = u s = or w 2 = 0. Then eq. (21) yields 
V(p — p s ) = and 



'dP\ 



\dT J Fp v-v. 
where s s = s s (T, P, 0) and 



S ^ : (24) 



d» B (T, P, 0) _ tT dr ds s (T',P,0) _ (25) 



dP Jt x dP 

The quantity v s is rather small and may be neglected in (24). From (14) 
we obtain T x (ds a /dP) = —T(ds a /dT) dT x /dP. Using this, dT x /dP w 
— O.OlK/bar and the entropy s s of Fig. 1, the integral (25) can be evalu- 
ated numerically. The resulting \v s /v\ has a similar temperature dependence 
as s s /s. The absolute values of \v s /v\ are much smaller, 

< 3 • 10" 4 . (26) 

The term v s has not been considered in previous generalizations [O, ITTJ of 



(19). It is derived theoretically; it follows from curlu s = and the pressure 
dependence of the superfluid entropy. 

2.3 Summary 

The two-fluid equations with s s ^ are 

d t p + V(p n u n + p s u s ) = 0, (27) 

d t (ps) + V(p n s n u n + p s s s u s ) = 0, (28) 

d t (p n u n + p s u s )i + dj (P5ij + p n u ni u nj + p s u si u Bj ) = , (29) 

md t u s + m(u s V)u s = -V(p - p s ) (30) 

where 

p s (T, P, w 2 ) = - [ T dT' s a (T', P, w 2 ) . (31) 

For s s = these equations reduce to (1) - (4). They are still eight equations 
for eight variables. The superfluid entropy s s is just a further macroscopic 
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quantity; as any other Galilean invariant macroscopic quantity it is of the 
form (5). It does not constitute a new independent variable. 
In addition we note: 

1. The equations (27) - (30) have to be supplemented by terms describ- 
ing dissipative effects. These terms might also contain corrections of 
the size 0(s s /s) relative to their well-known form ||. For practical 
purposes (like an estimate of the damping of sound modes) we will use 
unmodified dissipative terms. 

2. Eq. (6) follows from (30). Onsager's quantization rule is unchanged. 

Together with the two-fluid equations the underlying microscopic conception 
is slightly modified (section 1.2). As usual there is a macroscopic wave func- 
tion if) = ^/po" exp(i$(r)) which determines the supervelocity u s = ?iV$/m. 
What is given up is the identification of \ip\ 2 with p s . At the same time Lon- 
don's main point, namely the relation between the IBG and liquid helium, is 
reinforced by reconciling (7) with (8). 



3 Sound modes 
3.1 Introduction 

A major and important application of hydrodynamic equations is the eval- 
uation of sound modes. This application is of particular interest for testing 
the modified two-fluid model because sound velocities can be measured with 
high accuracy. We consider the first and second sound of bulk He II (section 
3.2) and the fourth sound of clamped He II (section 3.3). We determine the 
corrections in the sound velocities due to a non- vanishing superfluid entropy. 
The derivation of the fourth sound for s s 7^ leads nearly automatically to 
a new sound mode which we call sixth sound. The detailed calculations are 
given in Ref. ||14|| . The sixth sound has already been presented in a short 



letter [15 



The standard ansatz for sound modes 

T(r,t) = T + AT exp(i(kr — cut)) , (32) 

P(r,t) = P + AP exp(i(kr-u;i)) , (33) 

u n (r,t) = u ni0 + Au n exp(i(kr - ut)) , (34) 

u s (r,t) = u S;0 + Au s exp(i(kr - tut)) (35) 

is inserted into (27) - (30). The constant values T , Po, u nj0 and u Si0 solve the 
equations. Quadratic and higher order terms in the (small) amplitudes AT, 
AP, Au n and Au s are omitted. This yields a linear, homogeneous system 



9 



of equations for the amplitudes. For a non-trivial solution the determinant 
of the coefficient matrix must vanish. This condition yields solutions of the 
form uj v = uj v (k) (for a specific direction of k). The ratio 

C = ^ (36) 
k 

is the velocity of the sound wave. The sound velocities of the common two- 
fluid model (1) - (4) are denoted by c v $. Because of 



— < 2 • 10" 

s 



< 3 ■ 10" 4 (37) 

the differences between c v and c v $ are expected to be small. 

All coefficients in the linearized equations are taken at T , Pq and Wq = 
( u n,o — u s,o) 2 - The final results are expressed by thermodynamic quantities 
at wo = 0. Eventually we omit the index zero and use the notation 

X(T,P)=X(T o ,P o ,0). (38) 



3.2 First and second sound 

The equilibrium state of bulk He II has a given temperature T , pressure P 
and 

u n ,o = u Si o = . (39) 

We insert (32) - (35) into (27) - (30) and determine (36). This calculation 
is quite analogous to that 0, [| in the common two-fluid model. Therefore, 
we restrict ourselves to a presentation of the results. Because of (39) the 
coefficients in the linearized equations and the sound velocities are of the 
form (38). 

For the first and second sound the leading corrections to Ci j0 and c 2 are 
given by 



*<T,io=^(i-*(i-!*)4-*^4 m 

\ s V cpj u{ v 2pcp \oT J 1 



(40) 



and 

c 2 (T,P)=c 2 , (l--) • (Ml) 

We use the abbreviations 



^ and ^Jlh*-. (42) 
V op \Jmp n c v 
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By cy and cp we denote the specific heats at constant volume and pressure, 
respectively. The standard approximations for c li0 and c 2j0 are 



ci, ~ Mi and c 2 , ~ «2 J — • (43) 

V c p 

The corrections to these approximations are discussed in detail in Ref. ||]; 
they are of the relative order (u%/ul) (1 — cy/cp) ~ 10 -5 . 

There are two correction terms in (40). The first one has the relative size 
2 • 10~ 7 , the second one 2 • 10~ 8 . These corrections are of academic interest 
only; they are even smaller than the error in c^o ~ u%. 

The corrections for the second sound are of the order s s /s < 2%; they 
might be observable. Using (24) we may define a 'fountain pressure' entropy 
s FP by 



dT) 



(44) 

FP 



This yields a common expression for c 2 and c 2i o: 

ymp n c P \ C2 (s s ^0). 



(45) 



3.3 Fourth and sixth sound 

In this subsection we consider He II in which the normal phase is clamped, 

u n (r,t)=0. (46) 

This reduces the number of variables to five, for which we may choose T, P 
and u s . Their dynamics is determined by the five equations (27), (28) and 
(30); eq. (29) is effectively replaced by (46). 

A standard device for measuring sound modes in clamped helium is a 
ring filled with powder (Fig. 2). The position (middle line) of the ring may 
be described by 

r ring = (R cos (p, R sin 0, 0) , = 0, . . . , 2tt . (47) 

The thickness of the ring is assumed to be small compared to the radius R. 
Then the r-dependence of the considered fields reduces to a 0-dependence 
and the supervelocity is parallel to the ring, 

u B (r,t)=u B (</>,t)e+. (48) 
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Powder-filled ring with He II 




current 



Figure 2: A ring with clamped helium is a standard device for fourth sound 
and persistent current experiments. We consider the solutions of the lin- 
earized equations of motion for this configuration. For s s = these solutions 
are the fourth sound and a static fountain pressure (FP) gradient AP/AT. 
For s s 7^ (and in the presence of a persistent current) the FP gradient 
becomes a sound mode which we call sixth sound. 



This reduces the number of variables and equations from five to three. Using 
(46) and (48) the field equations (27), (28) and (30) become 

d t p+^^(p s u s ) = 0, (49) 
1 d 

dt(ps) + - — (p s s s u s ) = 0, (50) 

md ' u ' + RT4,\^ ll ' + ^r) = °- (51) 

The following derivation is simplified by using the variables 

o 

TflU 

T, p c = p- p s + — — - and j = p s u s (52) 



rather than T, P and u s . Taking into account (46), the specific geometry, 
and the new variables, the ansatz for sound modes reads 

T((p,t) = T + AT exp(i[kR(p — cot]) , (53) 

p c ((f),t) = p cfi + A/i c exp(i[/ci?0 - cot]) , (54) 

J'(<M) = Jo + Aj exp(i[kR(f) - cot]) . (55) 

The equilibrium state has the constant values T , p Ct o and 

Jo = Ps,ou s , = const. . (56) 
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A persistent current (j 7^ 0) is a metastable equilibrium. 

Because of T(0, t) = T(<f> + 2ir, t) the possible /c-values are restricted to 



h — — 
k ~ R 



where n G {±1, ±2, ±3, . . .}. 



(57) 



Inserting (53) - (55) into (49) - (51), the resulting linearized equations are: 



UJ 



dp 
df 



UJ 



dp 
dp c 



ds ds s 
up — -kj 



V 



dT 



uj u s 



dT 



dp s 
dT 



ds ds s 

k 3 ^~ 

dp c dp c 

uj u s h p s k 

dp s 



d P 1, 
dj 

ds ds 
UJp ~dj + ( S ~ Ss ^ k ~ k i ~dj 



UJ U s 



dp s 
dj 



UJ 



W Ar \ 



A/i, 
Aj 



/ 



V J ) 



(58) 

The first line represents eq. (49), the third line eq. (51). For the second line 
we subtracted eq. (49) times s from (50). All coefficients in (58) are to be 
taken at T , // Ci0 and j . Here and in the following we suppress the index 
zero. In each partial derivative two of the three variables T, /x c and j are 
kept constant. 

For a non-trivial solution the determinant of the coefficient matrix in (58) 
must vanish. This condition yields three eigenvalues uo u (for given k), or three 
sound velocities c v . The three sound velocities may be distinguished by their 
dependence on the velocity u s : 



c„ = 



±c 4 oc l + 0(u s ) 
c 6 oc u s (l + C(m s 2 )) • 



(59) 



The solutions ±04 are those of the well-known fourth sound. For s s 7^ and 
u s 7^ we obtain a new sound mode which we call sixth sound. 

In the following we restrict ourselves to the leading order in u s . Higher 
order contributions are necessarily small because |w s | must be less than the 
critical velocity -u cr i t . In a specific experiment, these contributions can be 
further suppressed by choosing \u s \ -u cri t- The coefficients appearing in the 
linearized equations are written as 



X(T, /i c ,j 2 )=X(T, Mc) + OK 2 ) 
The condition for a non-trivial solution of (58) yields 



(60) 



c 4 



Ps 

mp 



dp ds 
{s - Ss) dT +P dT 



dp ds 
~dp~ c df 



dp ds 
df dp~ c 



(61) 
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and 



Ps 

03 = u s — 



ds s 
dT 



p ds s — s s dp 
dT + p dT 



(62) 



Since we restrict ourselves to the leading order in u s all quantities are to be 
taken at T, p c and j = 0, eq. (60). We switch now to the arguments T, P 
and w 2 = by writing 



dX 

dX 
dp c 



dX(T,p c ) 
dT 

dX(T,p c ) 
dp c 




V — V. 



(63) 
(64) 



We use (dp c /dT)p = — (s — s B ) and (dp c /dP)T — v — v s which follow from 
(22) and (25). With (63) and (64) we obtain from (61) and (62): 



and 



\ p V w 



1 + 



T(s- 2 - t> s /t> / 9p x 
pc P 1 - v s /v \<9T / 



p J 



Mi 



Pn f x _ £s 



C 6 = M s 



g - g s / ds s \ 
dT J P v-v B \dPj T 



p fds\ 
\dTjp 



1 + 



PnUj (1 ~ S s /s) 

p s ii? (1 - v s /v) 



+ 



s s 2 - u s 



p 1 — u s /w \dT J p 



(65) 



(66) 



These results are exact to any order in s s /s and f s /f , and to leading order in 
u s . We present now approximate but simpler formulae. Neglecting the higher 
order corrections eq. (65) becomes 



a(T,P) 



\ p 



Ps _ V> 




p 



tin 



(67) 

For v s = and s s = this result is well-known 0. The first term under 
the square root is the dominant one; therefore it is appropriate to retain the 
small correction v s /v in this term (as compared to the larger correction s s /s 
in the second term). The correction in the second term corresponds to that 
in (41); the one in the first term is different from that in (40). 

The sound velocity Cq is proportional to s s . For a more handy expression 
we may therefore neglect all corrections to Cq of the order of one percent. We 
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go back to (62) and use dsJdT = (ds s /dT)^ c « {dsJdT)^ and ds/dT 
(ds/dT)^. Because 



s dp 
~p df 



10" 



ds 



dT 



(68) 



the second term in the denominator in (62) can be omitted. Thus we obtain 



where 



P c M 



= T | | and c,, s = T 



<9T 



(69) 



(70) 



are the specific heats at constant chemical potential. For the practical eval- 
uation (next section) one may use the specific heats at constant pressure 
instead. 



3.4 Discussion of the sixth sound 

The fourth sound is experimentally and theoretically well-established. The 
following discussion centers therefore on the sixth sound. We derive the am- 
plitudes of this sound mode and its damping. 

For s s = the frequency of sixth sound becomes zero; the solution (32), 
(33) represents a static temperature and pressure gradient. For the discussion 
of this limit we may go back to the equations of motion (49) - (51): For 
d t — the first equation is solved by d(p s u s )/d(f) = 0, and the second one is 
trivially fulfilled. The third one yields p — p s = const. + 0{u^) or p m const.. 
That means that the sixth sound reduces to a static fountain pressure (FP) 
gradient in the limit s s = 0. This situation is changed for s s ^ 0: The 
second term in (50) is now non-zero (provided that u s ^ 0) and leads via 
the first term to a time-dependence. The originally static fountain pressure 
gradient becomes an oscillating mode; effectively, the FP gradient is shifted 
with velocity Cq along the ring. The sixth sound may be conceived as an 
entropy transport along the ring caused by the persistent current (section 4). 

As discussed, the sixth sound amplitudes solve eq. (51) by p ~ const.. 
This means 

(—) « - • (71) 

\AP J 6th sound S 
The corresponding ratio for the fourth sound is well-known (for s s = 0). 
Neglecting corrections of the order 10~ 2 this ratio is 



APJ 



4th sound 



PnU 2 2 Ts_ ( dp_ 
p s ul pep \dT / 



p J 



V -« V -. (72) 

s s 
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The fourth and sixth sound can be distinguished by their amplitudes and 
by their velocities: The temperature amplitude of the fourth sound is much 
smaller than that of the sixth sound (for a given pressure amplitude). The 
velocity of the sixth sound has a characteristic proportionality to u s . 

In order to calculate the damping of the sound modes the two-fluid equa- 
tions are supplemented by (unmodified) dissipative terms ||. The detailed 
calculation |14| shows that only the heat conductivity k contributes to the 



damping of the sixth sound. We derive this damping in a simplified way. 
Taking into account the heat conduction the first coefficient in the second 
row of (58) becomes 

ds , ds s run , 9 . . 

ujp — -kj—?-+i — k 2 ttO. 73 
dT dT T K J 

For the sixth sound with its dominant AT-amplitude in (58) this coefficient 
must be approximately zero. For k = this condition yields the velocity 
cq = uj/k of (69). For k ^ we obtain instead 

uj = c 6 k — i = upp — i r F p • (74) 

This result is of the same accuracy as (69). The real quantities co>fp and Tpp 
will be used and discussed in section 4. 



3.5 Detectability of a superfluid entropy 

Each of the calculated sound velocities has some corrections due to s s ^ 0. We 
discuss whether a possible superfluid entropy may be detected by measuring 
the sound velocity: 

First sound: The relative difference between c\ and c^o is of the order 
10~ 7 . This is many orders below the experimental accuracy (and also 
below the theoretical accuracy of Ci j0 ). 

Second sound: The velocity c<i has been measured with a systematic 
error of less than 0.4% ||. Assuming the validity of (19) this measure- 
ment has been interpreted as that of the entropy. This entropy may be 
called 'fountain pressure' entropy sfp, eqs. (44) and (45). According to 
Fig. 1 we expect a 1 to 2% difference between sfp and the true (caloric) 
entropy s. Unfortunately the absolute value of s at the A-point is also 
uncertain by about 2%. 

Improving the present experimental accuracy by some factor (say 5) 
one should be able to see a difference between s and spp, in particular 
because of the steep rise of s s just below the A-point (Fig. 1). This 
possibility has been discussed in more detail in Ref. [[ll[ . 
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Fourth sound: The fourth sound velocity (67) is dominated by the 
first term under the square root. The correction in this term is of the 
order 10~ 4 ; it is smaller than the experimental accuracy. 

Sixth sound: The sixth sound velocity is proportional to the superfluid 
entropy. The sixth sound is therefore the prime candidate for detecting 
and measuring the superfluid entropy. Its observability will be discussed 
in detail in section 4. 

The attempt of measuring s — spp by the second sound will not become 
obsolete by a possible detection of the sixth sound. The reason is that the two 
modes are sensitive to different modifications of the two-fluid equations. The 
sixth sound velocity is due to the s s -contribution in the entropy continuity 
equation (28); the /vterm in (30) yields higher order corrections to c 6 only. 
On the other hand, a main correction in the second sound velocity stems 
from the p s -term in (30). 

4 Observability of the sixth sound 

We start this section with a simplified, alternative derivation of the sixth 
sound. We discuss then in detail the possible observation of the sixth sound. 

4.1 Entropy transport by the sixth sound 

We assume that the ring of Fig. 2 carries a persistent current with the velocity 
u s , and that at a certain time there is the following temperature and pressure 
variation along the ring: 

5T((f)) = A cos(n0) , 8P((f>) « ^ 5T , n e {1,2,. . .} . (75) 

This variation implies p, ~ const, (the small l.h.s. of (30) and the p s -term 
may be neglected for the present purpose). The variation (75) is a fountain 
pressure gradient; for s s = it is metastable. 

Eq. (75) implies a corresponding variation of the entropy density: 

P 5s=^ 5T{4>) . (76) 

Due to (68) we may use p ~ const.. The continuity equation (49) implies 
then p s u s const.. 

For s s 7^ the persistent current carries entropy. For constant T and P the 
net entropy current Sj s vanishes. For (75) we obtain (using p s u s ~ const.): 

Sj s = p s m s Ss s = p s u s -^- 5T{4>) . (77) 
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Due to this heat current the entropy variation (76) is shifted along the ring. 
This shift takes place with the velocity 

c 6 = — ^ = u s . (78) 

pos p c M 

In this rather simple way the sixth sound may be understood as an entropy 
transport phenomenon. Section 3 shows that it is also - in accordance with 
the common nomenclature - a sound mode. 

The net entropy current (77) must change the entropy variation of (76) 
of the whole ring (consisting of the container, the powder and He II). For 
the considered experiment, the entropy density ps in (50) must therefore 
include the entropy of the powder and of the container. In practice, this can 
be taken into account by equating c M in (78) with the specific heat of the 
whole ring, c M = c r i ng - At the considered low temperatures, the specific heat 
(per atom) of helium is at least four orders of magnitude larger than that of 
normal solid material. Therefore, we may set c ring ps c(He). In addition, the 
specific heat at constant chemical potential may be approximated by that at 
constant pressure (implying an error of the order 1CT 2 ): 

c M = c ring ps cp(He) , k = K ring ps re(He) . (79) 

These arguments apply similarily for the heat conductivity k. 



4.2 Shift of the temperature variation 

Any temperature and pressure variation in an actual ring experiment is likely 
to have a FP component. In the presence of a persistent current this variation 
is shifted with the velocity c 6 along the ring. This shift should be observable. 
It has characteristic properties: 

• The shift velocity is proportional to that of the persistent current. 

• The shift velocity is independent of the amplitude of the temperature 
variation (as long as the linear approximation works). 

The observation of the first property would leave no doubt that the persistent 
current transports entropy, implying that the superffuid entropy is non-zero. 
The experiment would then determine this superffuid entropy quantitatively. 
If no such shift is observed the experiment would yield an upper limit for S s 
which is considerably lower than the present one. 

For a potential experiment it is useful to know the expected size of the 
effect. We use c^ s ps cp iS which implies an error of the order 10" 2 . We insert 
this and (79) into (78): 

c 6 p s T ds s (T,P) 



u s p cp dT 



(80) 
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Figure 3: Velocity Cq of the sixth sound in the ring experiment of Figure 2 
as a function of the temperature T for saturated vapour pressure. 



We use the prediction s s (T,P) of Fig. 1 and the experimental p s /p and cp. 
Fig. 3 displays the resulting temperature dependence of c 6 /u s . 

The shift velocity of the temperature variation (75) along the ring 
amounts up to a few percent of the supervelocity. Depending on the tem- 
perature range the shift is parallel or antiparallel to the persistent current 
(Fig. 3). 

In an experiment, the temperature will be monitored at fixed points of 
the ring. At a given point the temperature amplitude oscillates with the 
frequency 

u F p = n . (81) 



For s s = (or u s = 0) the considered mode becomes a static FP gradient 
along the ring. Therefore, we call (81) the 'fountain pressure' frequency. 

As an example we insert the values u s = 2cm/s and R = 2 cm of an 
actual experiment [[llj. Using cq/u s ~ 3% and n = 1 we obtain an oscillation 
period 2it/ujfp of about three minutes. 
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Figure 4: Number n osc of observable fountain pressure oscillations for a ring 
experiment (R = 2 cm, w s = 2cm/s) as a function of the temperature at 
saturated vapour pressure. 



Due to the damping (74) of the sixth sound the amplitude of the FP 
oscillation will be reduced by a factor e after 



[LOppI 



2nT 



FP 



2irmK 



ds s (T,P ) 

or 



Ms 



R 



n 



(82) 



cycles of oscillation. We insert n = 1, R = 2 cm and u s = 2 cm/s. We use the 
prediction s s (T,P) of Fig. 1 and the experimental p s . The heat conductivity 
k = «;(He) of Hell is not well-known. For evaluating (82) we assume k = 
0.05 J/ (msK); this value is cited in [[J as an upper limit for T = 2.1 K. The 
result is displayed in Fig. 4. 

At the maximum of s s in Fig. 1 the velocity Cq and the number n osc vanish. 
The exact position of this point is subject to the uncertainty in the predicted 
s s . It is also slightly shifted by the approximation c MjS ~ cp jS . 

The condition n osc ^> 1 for easy observability is fulfilled at most tempera- 
tures. Moreover, the number n osc oc \u s \R of observable oscillations could be 
increased by using a higher supervelocity or a larger ring. The sixth sound 
or, equivalently, the FP oscillations should be readily detectable. 

The integration of eq. (80) yields 



T dT ,^± 
A Ps/P 



Cp_ 

V 



The quantities cq/u s , p s / p and cp are functions of T' and P. All these quan- 
tities can be measured. For the evaluation of the integral they have to be 
measured in the temperature range from T\ (where cq/u s and p s /p vanish) 
to T and at fixed pressure. In this way the proposed experiment determines 
the superfluid entropy s s (T, P). 
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If the sixth sound is not observed the considered experiment yields an 
upper limit for the superfluid entropy which is roughly given by 



S J upper limit ^-osc \ S J predicted ^osc 

In this way the present upper limit could be lowered by at least one order of 
magnitude. 



5 Concluding remarks 

A superfluid entropy of relative size S s /S ~ 1% is not excluded experimen- 
tally. We have presented a straight-forward generalization of the two-fluid 
equations for a non-vanishing superfluid entropy. The investigation of sound 
modes leads to a specific proposal by which a superfluid entropy of the con- 
sidered size could be detected. A negative experiment would yield a new 
upper limit for the entropy content of the superfluid fraction. 
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